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ABSTRACT 

We have previously found a temperature-dependent upflow in the dimming 
region following a coronal mass ejection (CME) observed by the Hinode EUV 
Imaging Spectrometer (EIS). In this paper, we reanalyzed the observations along 
with previous work on this event, and provided boundary conditions for model- 
ing. We found that the intensity in the dimming region dramatically drops within 
30 minutes from the flare onset, and the dimming region reaches the equilibrium 
stage after ~1 hour later. The temperature-dependent upflows were observed 
during the equilibrium stage by EIS. The cross sectional area of the fluxtube in 
the dimming region does not appear to expand significantly. From the observa- 
tional constraints, we reconstructed the temperature-dependent upflow by using 
a new method which considers the mass and momentum conservation law, and 
demonstrated the height variation of plasma conditions in the dimming region. 
We found that a super radial expansion of the cross sectional area is required 
to satisfy the mass conservation and momentum equations. There is a steep 
temperature and velocity gradient of around 7 Mm from the solar surface. This 
result may suggest that the strong heating occurred above 7 Mm from the solar 
surface in the dimming region. We also showed that the ionization equilibrium 
assumption in the dimming region is violated especially in the higher temperature 
range. 

Subject headings: MHD — plasmas — Sun: corona — Sun: flare — Sun: UV Radia- 
tion 
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INTRODUCTION 



A solar flare is often explained as a rapid energy conversion phenomenon from magnetic 
field energy to plasma energy. The released energy by a flare is so huge, and the total 
amount of energy often reaches to 10 32 erg within an hour. Over the past few decades, 
considerable effort has been devoted toward understanding the physical mechanism of solar 



flare, and various models have been proposed (e.g., Carmichael (1964); Sturrock (1966); 



Hirayama (1974); Kopp & Pneuman (1976)). Nowadays, it is widely accepted that magnetic 
reconnection is regarded as the fundamental energy conversion mechanism of eruptive flares. 
So far, various features expected from the model have been discussed and confirmed by the 



modern observations (e.g., Tsuneta et al. 1992 Ohyama & Shibata 1998 Yokoyama et al. 



2001 Teriaca et al. 2003). 



Another important aspect of the solar flare is its impact on surroundings. Most powerful 
flares are associated with a Coronal Mass Ejection (CME) (e.g., Kahler 1992 Yashiro et 



al. 2006). Because a CME releases a huge amount of plasma and magnetic field into 



the interplanetary space in a short duration of time, the equilibrium in the solar corona 
breaks down. Therefore, large-scale magnetic reconfiguration in the solar corona occurs. 
The opposite scenario also can work; a large-scale magnetic reconfiguration can cause a 
CME. Although the causal relationship between large-scale magnetic reconfigurations and 
CMEs is still not clear, a solar flare have a large impact on its surrounding through CME or 
large-scale magnetic reconfiguration. 



Coronal dimming in the EUV and/or soft-Xray range (e.g., Rust 1983 Sterling & 



Hudson 1997 Thompson et al. 1998), which is due to a loss of coronal plasma, is often 



observed as one of the main on-disk signatures of a CME. The dimming can persist longer 
than a day in the some extreme cases, and its mass loss contribution can reach more than 50% 



of the total mass of the CME (e.g., Sterling & Hudson 1997). There are two scenarios that 
can explain the origin of dimming regions. One is the eruption of a flux rope, and the other is 



interchange reconnection between closed and open magnetic field lines. Sterling fc Hudson 



(1997) found a pair of dimming regions which are located at either end of a preflare sigmoidal 



(S-shaped) structure. They concluded that their findings are consistent with the source of 
the CME being a flux rope that erupted, leaving behind the dimming regions. lAttrill et 



al. (2006) discussed the magnetic field topology changing associated with a CME. In this 



scenario, the magnetic loops of the CME expand and eventually push against the oppositely 
oriented open magnetic field of the polar coronal hole. This triggers successive magnetic 
reconnection, and some part of the magnetic loop topology changes from closed to open. 
Although it is not clear which is the dominant process, in either case the dimming region is 
likely to be a footpoint of open field line. The changing magnetic configuration from closed to 
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open within a short period causes the evacuation of coronal material towards interplanetary 
space through the rarefaction wave, and plasma cannot refill the open magnetic field lines 
for a long time. 

Quantitative measurements of dimming regions are essential to understand its contri- 
bution to CME formation or development. The EUV imaging spectrometer (EIS) on board 



Hinode (Kosugi et al. 2007) provides EUV spectra with higher spectral and spatial resolu- 



tion than ever before (Culhane et al. 2007). The first measurement of dimming regions by 



EIS was the flares and CMEs of 13 and 14 December 2006 in NOAA 10930. Ha rra et all 



(2007) studied the outflows in the dimming region of the 2006 December 14 events, and they 



found for the first time that the strong outflow exists at the loop footpoints in the dimming 



region which is located far away from the flare site. Imada et al. (2007) analyzed the flare of 



13 December 2006 and found a temperature dependence of the outflows; the hotter plasma 
shows the faster upflow. They claimed that the temperature-dependent upflows are the con- 
sequence of coupling between heating and bulk acceleration along the line-of-sight (height 



direction). Jin et al. (2009) have analyzed both of these December 2006 events and found 



that the velocity also correlates with the photospheric magnetic field. Harra et al. (2011) 



discussed the magnetic flux from the source of the CME by using the Doppler velocity dif- 
ference and the photospheric magnetic field measurement. They found that their estimation 
was consistent with the magnetic flux in the magnetic cloud observed at the interplanetary 
medium (1 AU) by in situ measurement. 

In this paper, we explore further the temperature-dependent upflow in the dimming 
region observed by Hinode/ElS. The aim of our paper is to reconstruct the temperature- 
dependent upflow by using a new method which uses the mass and momentum conservation 
law to determine the height variation of the plasma conditions in the dimming region. This 
paper is organized as follows. In the next section, we discuss the observational properties 
in dimming region of the 2006 December 13 event. Section 3 is devoted to the modeling of 
temperature-dependent upflow. We discuss the time-dependent ionization effect in Section 
4. Summary and discussion are given in §5. 



2. Observation of dimming region 

In this section, we discuss the observational properties in the dimming region on the 
2006 December 13 event. We have reanalyzed the event to obtain quantitative information 
to reconstruct the dimming region which is discussed in the next section. We also discuss 
the previous observational results of the same event. 
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2.1. Instrumentation 



The EIS instrument onboard Hinode is a high spectral resolution spectrometer aimed 
at studying dynamic phenomena in the corona with high spatial resolution and sensitivity. 
While making a large raster scan with 1 arcsec slit from 01:12 UT to 05:42 UT on 2006 
December 13, EIS observed an X3.2 flare. The flare occurred at 02:14 UT accompanied by 
the halo CME. The flare itself has been studied in detail (e.g., Kubo et al.||2007 Imada et al. 
2008 Asai et al.|2008 Minoshima et al.||2009[ ). EIS data from the raster were processed using 
the EIS team provided software (EIS_PREP), which corrects for the flat field, dark current, 
cosmic rays, and hot pixels. Each spectrum was fitted with a Gaussian profile, and the line 
width and Doppler velocity were determined. The slit tilt was corrected by eis_tilt_correction. 
For thermal reasons, there is an orbital variation of the line position causing an artificial 
Doppler shift of +/- 20 km s -1 which follows a sinusoidal behavior. This orbital variation of 



the line position was corrected using the method described by Imada et al. (2007), and the 



instrumental effects are reduced to below +/- 5 km s \ as a result. 

For studying the time evolution of dimming region, we made use of the Extreme ultra- 



violet Imaging Telescope (EIT) on SOHO (Delaboundiniere et al. 1995). We used the EIT 
images in the 195 A filter. The EIT data have a cadence of 12 min and a pixel size of 2.6 
arcsec pixels. The exposure times are ~12.6 sec during the event. We processed EIT data 
by using eit_prep for the calibration. 



2.2. SOHO/EIT observation 

An EIT base difference image using the 195 A filter is shown in Figure 1. A pre-flare 
image at 02:00 UT was subtracted from the post-flare image at 03:36 UT. The flare site is 
located around (350, -150), and there is a deep dark dimming region in the far east from the 
flare site. We highlighted the dimming region in which the intensity drops by more than 60% 
of its original intensity with a red contour. The time evolution of the dimming region can be 
clearly seen in Figure 2. Most loops or diffuse structures inside the red contour were blown 
off within the first 30 minutes, and the bright point-like structures, which may be related 
to foot points of coronal loops, loose their intensity gradually. Figure 3 shows the time 
evolution of the average intensity in dimming region (inside the red contour in Figure 1). 
We can clearly see that the intensity is rapidly decreasing and reaching the quasi-equilibrium 



stage. Similar features were observed in the event on the next day. (e.g., Harra et al. 2007). 



Long term evolution of this dimming region is discussed by Attrill et al. (2010). They found 
that concentrated downflows develop during the recovery phase of the dimmings and are 
also correlated with the same magnetic elements that were related to outflows during the 
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quasi-equilibrium stage. 



2.3. Hinode/EIS observation 



EIS observed the dimming region from 04:00 to 05:30 UT, and obtained the spectral 
information in the several emission lines. During this period, the dimming region already 
reached the quasi-equilibrium stage (see Figure 3). In Figure 4, we show the intensity, 
velocity, and line width map of He II (256.32 A; logT max . = 4.70), Fe XII (195.12 A; 
\ogT max = 6.11), and Fe XV (284.16 A; \ogT max = 6.30) obtained by EIS. The color scales 
in the velocity map range from -150 to 150 km s -1 . The color scales in the line width (Full 
Width at Half Maximum) map are normalized by the median of line width in whole field of 
view. The median values of line width in He II, Fe XII, and Fe XV are 0.102, 0.073, and 0.085 
A, respectively. Note that EIS scans from west (right) to east (left). The time stamps of slit 
scanning are also shown in the top of Figure 4. In the intensity map, the bright flare site can 
be observed in western part, and we can clearly see the dimming region in the far east from 
the flare site, which is framed by white dashed lines. The temperature-dependent upflow in 



the dimming region, which was discussed in Imada et al. (2007) and Jin et al. (2009), was 



also clearly seen in the Doppler velocity map in Figure 4. The strong outflows (~ 100 km 
sec -1 ) in higher temperature are located in the boundary of the dimming region. The line 
width in the dimming region also shows the temperature dependence. This is because the 
inhomogeneities of flows in the dimming region produce the spatial variations of Doppler 
velocities, which result in an apparent broadening of the integrated profile. Therefore, it is 
reasonable that the line widths in the dimming region also have temperature dependence. 



The relationship between flows and line broadening are discussed in Imada et al. (2008) and 



Dolla fe Zhukov (2011) 



To be confident of our results we need to understand any possible contamination from 
line blends. Generally, Fe XII and Fe XV lines are very strong in the corona, and we do 
not need to take into consideration line blending. On the other hand, in the He II line, 
there are strong line blends with the Si X which has a formation temperature of ~ 10 6 2 K 



(Young et al. 2007). However, the blend is not significant in the dimming region, because 



the line intensity from the coronal temperature plasma is very weak (Figure 4). Hence, we 
can neglect the line blend for He II in the dimming region. The line blend in the He II line 



is also discussed in Jin et al. (2009), and they also concluded that the contribution of the 



Si X line is more significant in the active regions than in the dimming areas. 



In the regions where the Doppler velocity is the highest in Fe XV, we fit the line profiles 
using a double gaussian for several ions. The line profiles are clearly separated into two 
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components - the fast upflow and the stationary component (not shown here, see Imada 



et al. (2007)). The fast upflow may correspond to the magnetic field whose topology has 



been changed from close to open by the flare activity. The other stationary component may 
correspond to the magnetic field whose topology was not changed by the flare activity. It 
is plausible that these two components are mixed within our spatial resolution, because the 
dimming region locates very far from the flare site. Figure 5 shows the relationship between 
the velocity (fast) and temperature in the dimming region. The line formation temperatures 
are defined by the temperature that yields its peak abundance with the CHIANTI code 
(Landi et al. 2006). This naturally means that we assumed the ionization equilibrium to 



derive the line formation temperature. The line blending effect does not have a significant 
effect on the velocity estimation in the result (see, Imada et al. 2007). The result of the 



relationship between the temperature and upflow velocity is also discussed in Jin et al. 



(2009). The velocities in our result is much faster than those in Jin et al. (2009), because 



we fitted the line profiles by the double gaussian. To discuss the relationship between the 
temperature and velocity quantitatively, we fit it by the following equation; 



T{v) = a 2 v 2 + a,\V + a , 



X 



where, T, v, and a 2 _o are the electron temperature (K), upflow velocity (m sec -1 ), and fit 
coefficients. The solid line in Figure 5 shows the fitting result, and the fit coefficients are 
also shown in Table 1. The correlation coefficient (r) is 0.95. 

The area of the EIS dimming region at different temperatures in the corona are also 



studied by Jin et al. (2009). They estimated the area with an intensity drop larger than 5% 



at the eruption phase compared with the pre-eruption phase. They found that the dimming 
area weakly depends on its temperature. From their result, the area increases by a factor of 
1.6 between 0.4 and 1.6 MK plasma. They also discussed the relationship between magnetic 
filed and upflow, and found that the fast upflows are located at relatively strong magnetic 
fields. 

There are two interpretations for the temperature-dependent upflow in the dimming 
region. One is that the temperature and velocity depend on the height from the solar sur- 
face (line of sight direction). The other is that the different loops in the dimming region 
have different temperatures and velocities in sub-resolution size (the multi-strand model). 
Both of the scenarios can explain the temperature-dependent upflows in the dimming re- 



gion. The recent observations also support both of the interpretations, (e.g., Jin et al. 2009 



Robbrecht & Wang 2010 Warren et al. 2011). However, the multi-strand model generally 



causes large broadening, distortion, and humps in line profile, (e.g., Ugarte-Urra h Warren 



2011). Therefore, the line profile may not separate into the two distinct components (fast 



and stationary), which is observed by EIS in the dimming region. Further, the area of the 
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dimming region also weakly depends on the temperature. The observations might support 
the former scenario, because generally the fluxtube expands radially in the solar corona. 
For those two reasons, we assume that the temperature and velocity depend on the height 
from the solar surface in the dimming region, although we do not deny the possibility of the 
multi-strand model. 



2.4. Summary of observation 

A summary of the observations is given below. 

• The intensity (Fe XII) in the dimming region dramatically drops within 30 minutes 
from the flare onset. 

• The dimming region reaches the equilibrium stage after ~1 hour after the flare onset. 

• EIS observed the temperature-dependent upflows in the dimming region. 

• The dimming region is in the equilibrium stage during the EIS observing period. 

• The temperature dependence of the upflow is characterized by Equation 1. 

• The area of the dimming region weakly depends on the temperature (the expansion 
factor is ~1.6 between 0.4 and 1.6 MK plasma). 

3. Modeling 

3.1. Method for Modeling 

In this section, we try to reconstruct the temperature-dependent upflow and determine 
the height dependence of plasma conditions in the dimming region by using the observational 
facts discussed in the previous section. Because the reconstruction will be based on our 
observations, our model can only apply in the range of 16 < v < 160 km s _1 (0.4 < T e < 2 
MK). We assumed the temperature-dependent upflow is a quasi-steady, because the flow 
was observed in the equilibrium stage of the dimming region. The fluxtube is vertical to the 
solar surface, and the cross sectional area is expanding with height from the solar surface. 
The fluxtube contains low (3 plasma, and we assumed that cross-field motions are negligible 
so that a one-dimensional description is adequate. Neglecting time derivatives in the hydro 
dynamic equations produces the quasi-steady flow model. The mass conservation equation, 
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equation of state in fully ionized gas, and momentum conservation equation in one fluid 
description are given by 

mnvA = const, (2) 

p = 2nk B T, (3) 

dv 1 dp / R s \ 2 

ds ran ds \s + R s J ' 

where n, v, A, p, and T are the electron number density, bulk velocity, cross sectional area, 
total pressure and temperature, respectively, m is the mean mass for solar abundances (1.257 
times of the hydrogen mass), k B is the Boltzmann constant, go is the solar surface gravity, 
R s is the sun radius, and s is the height from the solar surface. We neglected the effect of 
He or the other minor elements in Equation 3. 

In order to describe the flow precisely, we rewrite Equation 4 using Equation 2 and 3 as 
follows; 

^ 2k B T\ dv + 2k B dT _ ( R s \ 2 + 2k B T 1 dA 
rav J ds m ds \s + R S J m A ds 

The equation is almost the same equation which was used for the discussion of the solar wind 



(Parker 1958) or the discussion of the stationary siphon flows in coronal loops (e.g., Craig 



& McClymont 1986 Orlando et al. 1995). In a Parker-type wind, the flow is accelerated 
by a pressure gradient. The pressure drop is typically due to decreases in both temperature 
and density. In the flow described here, T is increasing with v (see Figure 5). Therefore, we 
anticipate that acceleration requires rapid divergence of the fluxtube area and that the flow 
is strongly heated. 

Equation 5 has a critical point at v — [2k B T ' /m\ 1 ^ 2 . Fortunately, the upflow velocity in 
the dimming region is less than [2k B T /m} 1 ^ 2 , for example [2k B T '/m] 1 / 2 = 113 km sec -1 at 1 
MK. We substitute T(v) (Equation 1) into Equation 5, and rewrite it in conservative form 
as follows; 

\ (v 2 - vl) - 2 -t (|«2 (v 2 - ^ 2 ) + ai (v- v ) + a log j) + ^ (T - T ) 

= ~9oR 2 (-^ + ^) + 2 ~t & Ti§ds', (6) 

where v , T , s , are the velocity, temperature, and height at the base of the coronal loop. 
We have chosen the footpoint velocity and height to be 16 km s -1 and 2000 km, respectively. 
To can be determined by Equation 1. We simplify each term in Equation 6 as follows; 
Sli — Si2 + Sl3 = —Sri + Sr2- The left-hand side in Equation 6 (Sli — Sl2 + Sls) can 
be determined when the velocity is determined, and the first term of right-hand side can 
be determined when the height from solar surface is determined. Only the second term 
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of right-hand side (Srz) is unknown. Therefore, we need to make an assumption of cross 
sectional area (A) to evaluate Equation 6. 



3.2. Radial expansion for fiuxtube 

The geometry of coronal fluxtubes is believed to fan out from the network boundaries 



in the photosphere to the corona, forming a canopy like structure (Gabriel 1976). This 
geometry leads to a varying cross sectional area of every fiuxtube along the height from the 
solar surface. We assumed the cross sectional area expands radially, and it can be written 
as follows; 

Ms) = Ao(P^-)\ (7) 



s + R. 

where \i is the parameter for expansion. The uniform, sub-radial expansion, radial expansion 
and super-radial expansion in cross section can be expressed by/i = 0,0</i<2,/i = 2 
and 2 < /x, respectively. In all these cases, Sr% can be written as follows; 



m J so s' + R s 
Further, we can evaluate the integration part in Sr 2 with 

T ^f m VTW, ds ' > r^TT*? > r »£ 7Tb? < (9) 

where Tmax is the maximum temperature in our range. 

The value of each term in the left-hand side in Equation 6 is shown as a function of 
velocity in Figure 6. We can clearly see that the total of left-hand side terms always has 
positive values. In Figure 7, we show Sri as a function of height by the solid line. Note 
that the horizontal axis in Figure 6 and 7 are different. Let us evaluate Sr 2 in the case of 
radial expansion (/i = 2). The maximum and minimum value obtained from Equation 9 are 
shown in Figure 7 by the dashed and dotted lines. Thus, Sr 2 should be located between 
the dashed and dotted lines. This result indicates that the right-hand side (—Sri + Sr 2 ) is 
always negative. Therefore, there is no solution of Equation 6 in the case of radial expansion, 
because the left-hand side in Equation 6 is always positive. 
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3.2.1. Super-radial expansion: Case I (fi = 45) 

From the slope of Sri and Sf® at s — sq, we can estimate the lower limit of \i from the 
condition being the positive value in the right-hand side of Equation 6 (dS^/ds > dSm/ds)] 

a > —— yo s - ~ 44. 10 
^ " 2k B T s + R s 1 ; 

Thus we assumed the super-radial expansion (// = 45) for the cross sectional area (Case 
I). The dotted-dashed and dotted-dotted-dashed lines in Figure 7 show the maximum and 
minimum value of Sr2 (a* = 45) evaluating from Equation 9. This result indicates that the 
right-hand side is always positive. Therefore we can solve Equation 6 with \i = 45 and obtain 
the height dependence of plasma conditions in the dimming region. We numerically perform 
integration of Sr®. 

The results for the height dependences of various quantities in the dimming region are 
shown in Figure 8. The horizontal axis shows the height from the solar surface, and the 
vertical axes show the velocity, temperature, cross sectional area, and density. The cross 
sectional area and density are normalized by those at s = 2000 Mm {=sq). The upflow is 
almost constant until a height of 6 Mm and dramatically accelerated up to 160 km sec -1 
from 6 to 9 Mm. The temperature also shows the same profile as the velocity because of 
Equation 1. The cross sectional area expands almost proportional to the height from the 
solar surface, and it reaches around 1.6 at 9 Mm. The density is gradually decreasing with 
the height, and rapidly decreasing from 6 Mm. The density ratios (n/no) at 7.5 and 9 Mm 
are almost 0.14 and 0.07, respectively. Note that the temperatures at 7.5 and 9 Mm are 
around 1.4 and 2 MK, respectively. 



3.3. Empirical Model for fluxtube 



The empirical modeling of the relative size of the fluxtube cross sectional area in the quiet 
sun was done by Chae et al. ( 1998[ ). They obtained the relationship between temperature 
and velocity from SOHO/SUMER observations in the quiet sun, and derive the estimation 
of the relative size of the fluxtube as a function of temperature. They also found the derived 



relative size of the fluxtube can be fitted to a functional form suggested by Dowdy et al. 



(1987 


) and 


Rabin 


(1991 



A{T) = A(T h 



1 
f 



(r 2 - 1) 



T 
Th 



'ID 
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with values of = 10 6 K, F = 31, and v = 3.6. We applied Equation 11 to Equation 6, and 
Sr2 can be written as follows; 

s »=£¥jC T "( i+(ia - i) ©y v <i2) 

Their discussion is concentrated on the down flows (< 10 km s _1 ) observed in the quiet sun 
transition region (below 1MK). Although their situation may be different from the dimming 
region, their results might include the characteristics of the cross sectional area. 



3.3.1. Strong dependence on temperature: Case II (v = 3.6) 



First, we have applied the parameters suggested by Chae et al. (1998) {T h = 10 6 K, T = 
31, and v = 3.6) to Equation 12, and solved Equation 6. The results of height dependence in 
the dimming region are shown in Figure 9. The figure format is the same as Figure 8. The 
upflow is linearly accelerated until a height of 150 Mm, and the velocity gradient becomes 
much steeper beyond 150 Mm. Roughly speaking, the temperature is increasing linearly in 
the entire region. The 2 MK plasma is located at 170 Mm. The cross sectional area expands 
almost proportional to the height from the solar surface, and it reaches around 25 at 150 
Mm. Because the velocity and cross sectional area are almost proportional to the height, 
the density is decreasing with height (~ s~ 2 ). The density ratios at 100 and 170 Mm are 
~ 0.015 and 0.004, respectively. It is also remarkable that the temperature at 100 and 170 
Mm are ~ 1.4 and 2 MK, respectively. 



3.3.2. Weak dependence on temperature: Case III (v = 1.5) 

The increase of cross sectional area between 0.4 and 1.6 MK plasma should be around 
~ 1.6, which we already mentioned in §2. It seems that the cross sectional area expands too 
much in the previous case (Case II). The above parameters are derived from the quiet sun 
observations. Therefore those may be different when we discuss the cross sectional area in 
the dimming region. Next, we assumed the parameters as Th = 10 6 K, T = 31, v = 1.5 for 
the relatively gradual expansion in the cross sectional area (Case III). The results of height 
dependence in the dimming region are shown in Figure 10. The figure format is the same 
as Figure 8. The characteristics in Figure 10 are very similar to that in Figure 9. Roughly 
speaking, the velocity, temperature, and cross sectional area are increasing linearly, and the 
density is decreasing with height (~ s~ 2 ). The difference between Figure 9 and 10 is the 
absolute value of the height from the solar surface. The maximum height is ~ 60 Mm in 
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Figure 10, although that is 180 Mm in Figure 9. The temperatures at 40 and 55 Mm are 1.4 
and 2 MK, respectively. The density ratios at 40 and 55 Mm are 0.063 and 0.03, respectively. 
Note that our result is not sensitive to Th or T but sensitive to v largely. We have tested our 
method with other parameters (not shown here). 

One may think that the increase of cross sectional area between 0.4 and 1.6 MK plasma 
is still large when we assume v = 1.5. Assuming a much smaller value in v (~ 1), the 
increase will become much smaller, and the height from the solar surface becomes much 
lower. However, when we assume v ~ 1, there are double valued solutions in the higher 
temperature range. Thus we cannot obtain the realistic solution with v ~ 1. 



3.4. Comparison between Cases I— III 



Let us discuss which is the most realistic result within our three examples. One of the 
biggest differences between the three results is the height from the solar surface. In Case 
I, 1.4 MK plasma (Fe XII) is located around 7.5 Mm from the solar surface. On the other 
hand, 1.4 MK plasma is located around 100/40 Mm from the solar surface in Case II /III, 
respectively. Unfortunately, we do not have any observational information of the height in 
this event. However, some similar events which occurred in the limb showed that the line 



formation height of Fe XII (1.4 MK plasma) is <10 Mm. Balasubramaniam et al. (2010) 
discussed the flare event on 2006 December 6. The flare occurred in the same active region 
as our event (NOAA 10930), and their characteristics are also very similar to our event 
(e.g., eruptive flare, GOES class, CME associated flare). They also found the upflow at the 
boundary of active region and dark region (transient coronal hole) in Fe XII (see Figure 16 in 
Balasubramaniam et al. (2010)). The upflow was observed only in the base of a fiuxtube (< 



10 Mm from the surface). These findings support our results in Case I. Further the increase 



of the cross sectional area also support Case I. The observation in Jin et al. (2009) indicated 



that the increases of cross sectional area between 0.4 and 1.6 MK plasma seems to be 1.6. 
Therefore, we conclude that the results in Case I is the most realistic. 



4. Time-dependent Ionization 

We have discussed the height variation of temperature-dependent upflow in the dimming 
region by using the steady hydrodynamic equation. To derive the relationship between 
the temperature and velocity from the observations, we assumed ionization equilibrium. 
Generally fast flow and rapid heating may cause non-equilibrium ionization, and so far many 
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papers discussed the time-dependent ionization with the models (e.g. 



Reale & Orlando 2008 Ko et al. 2010 Imada et al. 2011) and observations 



Dupree et al. 


1979 


ions (e.g., 


Kato et al. 



1998 Imada et al. 2009 Murphy et al. 2011). It is plausible that the ionization equilibrium 
assumption is violated in our situation, because both the temperature and velocity rapidly 
increase with height. In this section we try to evaluate the ionization equilibrium assumption 
in our result. 

In order to study the effect of transient ionization, we have calculated the time evolution 
of ion charge states. In our EIS observations, most of emission lines are from iron. Therefore, 
we concentrated on the time-dependent ionization of iron in this paper. The method we use 



to calculate the time-dependent ionization is the same as that described in Imada et al. 



(2011). The continuity equations for iron are expressed as follows; 
dnf e 



dt 



+ V 



Fe 

nf V 



n, 



Fe Fe 
n i+l a i+l 



(13) 



where nf e is the number density of the zth charge state of the iron, af e represents the 
collisional and dielectronic recombination coefficients, and Sf e represents the collisional ion- 
ization coefficients. The ionization and recombination rates were calculated using Arnaud 
& Rothenflug I (119851), lArnaud & Raymond I 019921), and iMazzotta et al. I (119981). Here 



we assumed that all ions and electrons have the same flow speed and temperature in the 
same position. The ionization and recombination coefficients (a and S) strongly depend 
on electron temperature and weakly depend on density. The timescale for ionization and 
recombination is proportional to n" 1 . 

To solve the time-dependent ionization, we need the absolute value of density. We set 
the density at the base of fluxtube (no) to 10 10 and 10 11 cm -3 . These values are reasonable 
for the density in the upper chromosphere or lower transition region (e.g 



Vernazza et al. 



3 ), and the time-dependent ionization calculation with the higher density 
= 10 11 cm" 3 ) are shown. The peak abundances of each charge state with 



1981). Figure 11 shows the height variation of the iron charge states in the super radial 
expansion model (Case I; Figure 8). From the top, the result of the ionization equilibrium 
calculation, the time-dependent ionization calculation with the lower density assumption 
(no = 10 10 cm" 
assumption (no 

the ionization equilibrium assumption are shown with vertical dotted lines. From Fe VIII to 
Fe X, the ionization equilibrium assumption seems to be reasonable. On the other hand, it 
is clearly seen that the ionization equilibrium assumption is violated from Fe XI to Fe XV 
in the case of lower density assumption. This is because there is a steep temperature and 
velocity gradient around 7 Mm. Even in the case of higher density assumption, the ionization 
equilibrium assumption is violated in Fe XIV and Fe XV. 



We also perform the same analysis in the empirical rapid expansion model (Case II; Fig- 
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ure 9). Figure 12 shows the height variation of the iron charge states in Case II. The figure 
format is the same as Figure 11. It is plausible that the ionization reaches the equilibrium 
stage in this mode, because the gradient of temperature and velocity in space is relatively 
gradual compared with the previous model. However, even in this case, the ionization equi- 
librium assumption is violated in Fe XIV and Fe XV. This is because the rapid expansion 
causes the low density in the high temperature range, and the timescale of ionization in- 
creases. Therefore, in either case, the ionization equilibrium assumption may be violated in 
the higher temperature range. 

The violation of ionization equilibrium naturally causes the increasing of the line for- 
mation temperature. For example, the peak abundances of Fe XII in the case of ionization 
equilibrium (top of Figure 11) is located around 7.5 Mm (yellow dotted line). Thus, Fe XII 
represents ~1.4 MK plasma, because the temperature at 7.5 Mm is 1.4 MK in Figure 8. 
However, the peak abundances of Fe XII in the case of time-dependent ionization (middle of 
Figure 11) is located around 8 Mm. This indicates that Fe XII represents ~ 1.7 MK plasma, 
because the temperature at 8 Mm is 1.7 MK in Figure 8. In the higher temperature range, 
the line formation temperatures systematically increase by considering time-dependent ion- 
ization. Therefore, some modification may be needed to the relationship between the velocity 
and temperature in Figure 5, especially in the higher temperature range. The relationship 
may become close to linear after considering the time-dependent ionization. 



5. Summary and Discussion 

We have discussed the temperature-dependent upflow in the dimming region observed 
by Hinode/ElS. We have reanalyzed the dimming event on 2006 December 13, and described 
the previous results for this event. The observations show; 1) the intensity in the dimming 
region dramatically drops within 30 minutes from the flare onset, and the dimming region 
reaches the equilibrium stage after ~1 hour later, 2) the temperature-dependent upflows are 
observed in the dimming region during the equilibrium stage, 3) the area of the dimming 
region also weakly depends on the temperature (the expansion factor is ~ 1.6 between 0.4 and 
1.6 MK plasma). We carried out the reconstruction of the temperature-dependent upflow 
by using our new method with the assumption that both the temperature and velocity 
depend on the height from the solar surface. The height variations of plasma conditions 
in the dimming region are demonstrated for three different cases, and we discussed which 
result is the most realistic. What we found from the reconstruction is as follows; 1) the 
super radial expansion of the cross sectional area is required to satisfy the mass conservation 
and momentum equations, 2) the temperature and velocity gradient with height become 
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steep around 7 Mm from the solar surface, 3) the ionization equilibrium assumption may be 
violated in the higher temperature range. What we found is summarized in Figure 13. 

In our analysis, we assumed that both the temperature and velocity depend on the 
height from the solar surface. We already mentioned that the other scenario (the multi- 
strand model) can also explain the observation. Note that it is still consistent even if we 
apply our result to the multi-strand model, because the higher temperature loop (~ 2MK) 
also contains the cold component (~ 1MK) at the base of the loop. The difference between 
two scenarios is that the heating and acceleration are stopped before reaching to 2MK in 
some loops. 

Let us discuss the impact of temperature gradient to the flow acceleration. To under- 
stand intuitively, we assume p oc e~ s and T = \±s + A2, where Ai t 2 is the positive parameters. 
After substituting p oc e~ s and T to Equation 3, we can derive n oc (Ais + A2) _1 e~ s . Then we 
found that the first term in right-hand side in Equation 4 ((— 1 /mn)dp / ds) is proportional to 
Ais + A 2 . Roughly speaking, the second term in right-hand side in Equation 4 is negligible. 
Therefore, the steep temperature gradient causes the rapid acceleration. The steep temper- 
ature gradient causes the steep density drop. This is the reason why the strong temperature 
gradient causes the strong acceleration. Figure 14 shows each term in Equation 4 in Case I, 
and we can find that the trend of (—l/mn)dp/ds is very similar to that of the temperature 
in Figure 8. 

We also perform the time-dependent ionization calculation in our models, and found 
that the ionization equilibrium assumption may be violated in the higher temperature range. 
This naturally causes the increasing of the line formation temperature, and we need some 
modification to the relationship between the velocity and temperature in Figure 5. Further, 
the modification causes some another modification in our reconstruction of the dimming 
region. It may be possible to find the reconstruction solution which can satisfy both the 
hydrodynamic equation and the time-dependent ionization. However the time-dependent 
ionization process heavily depends on the absolute value of density. Unfortunately, we do 
not have any pairs of emission lines for density diagnostics in our observation. Thus we 
cannot conclude that the modification is really needed or not at this stage. 

We now discuss the energy balance in our model. The energy equation in quasi-steady 
flow can be written as follows; 

3 1 dT dv 1 dA 1 d ( ,*dT\ 9 , . , tt , <N 

where kq, A(T), and H are the thermal conductivity along the magnetic field line, the 
radiative loss function, and the heating function. From our reconstruction of the flow in 
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the dimming region, we can derive most of the terms in Equation 14. The exceptions are 
the absolute value of density and the heating function. Thus we can estimate the heating 
function with the assumption of density at the base of fiuxtube. We have calculated the 
heating function with no = 2.5 x 10 10 cm -3 in Case I. Figure 15 shows the energy balance 
in the dimming region. The solid (black), dashed (blue), dotted (red), and dotted-dashed 
(green) lines show the height variation of the heating function, radiative cooling, thermal 
conduction, and left-hand side of Equation 14, respectively. At the bottom of the fiuxtube, 
the radiative cooling is dominated because of the low-temperature and high-density plasma. 
Therefore the heating function should be large enough to maintain energy balance. There is 
a steep gradient in the velocity and temperature between 6 and 7 Mm in Figure 8. Thus the 
enthalpy transport and thermal conduction are dominated in this region. Above the 7 Mm, 
the temperature gradient becomes gradual. This naturally causes the rapid reduction of the 
thermal conduction, although the enthalpy transport reduces gradually. Thus an increase 
in the heating function is required. Our result in Case I suggests that the strong heating 
function may be located above 7 Mm. Note that most of Equation 14 is dependent on the 
absolute value of density, and the dependence on the density is different in the each term; 
radiative cooling oc n 2 , left-hand side oc n 1 , thermal conduction oc n°. Unfortunately we do 
not have the observational information of the density. Thus, we do not discuss the absolute 
value of heating function at this stage. 

Density diagnostics are crucial for the modeling of the dimming region. Consistent 
modeling, taking care of the time-dependent ionization and energy equation, of the dimming 
region with the density diagnostics is for future work. Another important question is the 
stability of the dimming flow in the category of hydro and/or magnetohydro dynamics. So 



far, the stability of coronal loops in many kind of instability were tested (e.g., Parker 1953 



Hood & Priest 1980 Zweibel 1981). The dimming flow seems to persist for a long period. 
The steady solutions dramatically change in some cases when the conditions are slightly 
modified. Further, the critical point seems very close to the solar surface, although it is not 
clear whether the flow pass through a critical point or not. It is physically important whether 
the flows become super sonic or not. This may also affect the stability of the dimming region. 
The stability analysis of dimming region is also future work. 
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Table 1: Fitting results. 

~q 2 (Ks 2 m- 2 ) ^(Ksm- 1 ) a (K) 
-6.00xl0~ 5 2.27X10 1 -2.05xl0 4 
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Fig. 1. — EIT base difference image with preffare image at December 13, 02:00 UT sub- 
tracted from the image on December 13, 03:36 UT. The red contour represents the dimming 
region. 
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Fig. 2. — Time series of the dimming region. The right figures are the magnified figure of 
the dimming region. The red arrows show the bright region inside the dimming region. The 
color scale was changed from the left figures. 




Fig. 3. — Time variation of the intensity in the dimming region (inside the red contour in 
Figure 1). The EIS observing time of the dimming region is also shown. 
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Fig. 4. — Intensity, velocity, and line width map of the dimming region observed by Hin- 
ode/EIS. The color scales in the velocity map range from -150 to 150 km s _1 . The color 
scales in the line width (Full Width at Half Maximum) map are normalized by the median 
of line width in whole field of view. The median line widths of He II, Fe XII, and Fe XV are 
0.102, 0.073, and 0.085 A, respectively. 
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Fig. 5.— Relationship between temperature and velocity in the dimming region. The 
velocities were estimated by Doppler velocities, and the line formation temperatures were 
derived from ionization equilibrium assumption. 
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Fig. 6. — Each terms in the left-hand side as a function of velocity in Equation 6. Solid, 
dotted, dashed, and dash-dotted lines show Sli — Sl2 + Slz, Sli, Sli, and Sl3, respectively. 
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Fig. 7. — Each terms in the right-hand side as a function of height from the solar surface in 
Equation 6. 
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Fig. 8. — Height variation of plasma conditions in the case of /i — 45 (Case I: super radial 
expansion). From the top to bottom, velocity, temperature, increase of cross section, and 
density drop are shown. The cross sectional area and density are normalized by that at 
s = s - 
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Fig. 9. — Height variation of plasma condition in the case of v — 3.6 (Case II: empirical 
strong expansion model). The figure format is the same as Figure 8. 
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Fig. 10. — Height variation of plasma condition in the case of v — 1.5 (Case III: empirical 
weak expansion model). The figure format is the same as Figure 8. 
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Fig. 11. — Time-dependent ionization in the case of fi = 45 (Case I: super radial expansion). 
From the top to bottom, the results with ionization equilibrium, non-equilibrium (no = 10 10 
cm -3 ), non-equilibrium (n Q = 10 11 cm -3 ) are shown, respectively. The vertical dotted lines 
show the heights of peak abundance in each elements in the ionization equilibrium results 
(Top). 
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Fig. 12. — Time-dependent ionization in the case of v — 3.6 (Case II: empirical strong 
expansion model). The figure format is the same as Figure 11. 
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Fig. 13. — Schematic illustration of our observation and reconstruction in and around the 
dimming region. T e shows the electron temperature in the dimming region. A and N e show 
the cross sectional area and electron density which are normalized at the base (2 Mm). 
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Fig. 14. — Each terms in Equation 4 (Case I: super radial expansion). 
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Fig. 15. — Energy balance in the dimming region (Case I: super radial expansion). 



